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^ ; Abstract 

"^ I A certain truncation of the overlap (domain wall fermions) contains k flavors of Wilson- 

Dirac fermions. We show that for sufficiently weak lattice gauge fields the effective mass 
of the lightest Dirac particle is exponentially suppressed in k. This suppression is seen to 
disappear when lattice topology is non-trivial. We check explicitly that the suppression 
holds to one loop in perturbation theory. We also provide a new expression for the free 
fermion propagator with an arbitrary additional mass term. 
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1. Introduction 

The overlap exactly preserves the non-anomalous global chiral symmetries of vector- 
like gauge theories [1] on the lattice. Relatively to chiral gauge theories the overlap sim- 
plifies considerably in the vector-like situation [2] but does not provide yet a practical 
alternative to more traditional simulations of QCD. The fundamental reason for the sim- 
plification is that in the vector-like case the Dirac operator can be viewed as a square 
matrix even in topologically nontrivial backgrounds. For chiral theories the shape of the 
matrix can be square or rectangular depending on the gauge background and therefore 
cannot be kept fixed: one needs an infinite number of fermions. For vector like theories 
the rectangular shapes of the two Weyl components always complement each other and can 
be combined in one rigid square shape. As a result, the overlap admits a truncation to a 
finite number of fermion fields (k). The truncated overlap can be incorporated in practical 
simulations of QCD and holds the potential to become competitive in the realm of small 
quark masses [3,4,5]. (The parameter k can also be viewed as the length of an extra dimen- 
sion connecting two "domain walls" on which the Weyl components live [6,7,8].) When k is 
taken to infinity one regains the vector-like version of the overlap where chiral symmetries 
are exact. At tree level, the truncated overlap has an exponentially small (in k) quark 
mass. As stressed recently in [3], the mechanism for this suppression is quite generic. The 
generic features indicate that the suppression ought to hold also after radiative corrections 
are taken into account. 

But, the indication does not constitute proof even by normal physics standards: It is 
worrisome that exact masslessness is protected at infinite k by an analytical index of a mass 
matrix, A4 [6]. The index, di'm{ker{M.^)) — dim{ker{M.)) is present at infinite k (infinite 
number of fiavors) and is robust under small perturbations as any index would be. But, 
for any finite /c, dim{ker[jvi'^)) = dim{ker[M.)) and there is no index. Thus a non-smooth 
behavior as /c ^ oo cannot be easily excluded. Chiral symmetries are notoriously difficult 
to maintain on the lattice and we must be careful. 

If any of the good features of the overlap are to be even partially preserved by the 
truncation we need new alternative ways to understand the suppression of the bare quark 
masses in the vector-like case, without relying on continuity in k at infinity. This prompted 
a search [3] for other interpretations of the mechanism: As was shown there, the mass 
suppression could be viewed as the outcome of a see-saw mechanism [9] or, alternatively, 
of an approximate conservation law associated with a chiral global symmetry in a way 
first devised by Froggatt and Nielsen [10]. In [3] it was also found that a convenient link 
between finite and infinite k is provided by a class of orthogonal polynomials associated 



with a certain measure on the real hne. The effect of the truncation is simply to keep only 
the first k polynomials, while the measure is /c-independent. Therefore, the single place 
k enters is in the index of the last polynomial kept. The coefficients of the polynomials 
themselves are ^-independent. This makes it possible to trace the complete dependence on 
k. Previous approaches [8,11,12] resorted to approximations in which parts of the objects 
needed at finite k were replaced by their infinite k limits. 

In the particular model related directly to the overlap version employed in [2] the 
polynomials turned out to be the Chebyshev polynomials of second kind. Other mass 
distributions would lead to other sets of orthogonal polynomials but would work similarly. 

Let us summarize the situation: None of the arguments in favor of an exponential 
suppression at finite k are fully compelling. The analytical index disappears, and even 
continuity in k is not a trivial matter. While the orthogonal polynomials analysis does 
indicate continuity, the exponential character of the suppression is evident only at tree 
level. The see-saw point of view more or less ignores doublers while the F-N picture relies 
heavily on chiral symmetries, a notoriously slippery concept on the lattice. Therefore, a 
direct check seems to be needed. 

As alluded above, we are not the first to undertake the task of checking approximate 
masslessness in a direct way [11,12]. Actually, an even earlier calculation to one loop order 
for the exactly massless case can be found in [13]. There, the perturbative calculation was 
carried out for the overlap in the more general, chiral, case. This calculation essentially 
checked whether the index argument survived one loop radiative corrections and concluded 
that it did. 

The work first announced in [11] and more fully discussed in [12] is still not complete, 
but it already claims to have established that masslessness is protected to one loop order. 
These authors ignore the relevance of the index at infinite k, and this blurs one of the 
main effects of the truncation. It is troubling that the calculation in [12] does not use the 
exact free fermion propagators but replaces them by propagators at /c = oo. The latter 
propagators have one strictly massless quark and therefore extra infrared singularities 
appear. Moreover, the k = oo propagator "knows" about the analytic index at A; = oo and 
it is exactly the disappearance of the index at finite k that makes the perturbative test so 
relevant. In short, the replacement of the exact propagators by the k = oo ones implies 
that certain /c-dependent terms have been ignored and it is not proven that the ignored 
terms are sub-dominant at large k. 

A full computation is somewhat tedious. We choose to avoid unnecessary details 
and focus on what the systematics are. We wish to know which features are essential for 
mass suppression and which are not. At tree level, the mechanism of suppression appears 



to be quite generic, and we would like to know how much extra "tweaking" one could 
do to improve it even further, without radiative corrections acting destructively. If the 
suppression indeed holds beyond tree level, it should do so for all kinds of variants of the 
gauge action or the fermion-gauge coupling. At tree level, even strict gauge invariance 
does not appear to be essential for mass suppression. Therefore we expect the suppression 
to work diagram by diagram (as long as we sum over all flavors in each - so the diagram 
"knows" about the full flavor structure) and even before the momentum integration is 
done, i.e. directly at the level of the Feynman integrand. 

Naively, the suppression appears quite miraculous. Consider flrst the basic example of 
a single Wilson-Dirac fermion on the lattice: In addition to the "light" Dirac fermion there 
are "heavy" ones, the doublers. The gauge field can turn a light right handed Dirac fermion 
into a heavy right handed doubler. The latter strongly mixes with its left handed heavy 
partner which, again via the gauge interaction, can become a light left handed fermion. 
Thus, the light Dirac particle gets a mass term of order g"^, where g is the gauge coupling. 
There also is more direct source of mass in the gauge sea-gull coupling to the light fermion 
via the Wilson mass term. 

Before considering several flavors let us ask why radiative mass generation is unavoid- 
able. Probably the most convincing form of the answer is the following: If the right and 
left handed fermions did not mix we would have exact U{l)v x U{1)a and we know this 
cannot hold for all gauge flelds because of instantons. Once we allow the helicities to mix, 
it appears impossible to assure masslessness except by flne tuning (i.e. by adding a linearly 
divergent mass counter-term - the divergence is 0{a~^), where a is the lattice spacing). 

The truncated overlap contains a larger number (16A; — 1) of heavy Dirac fermions. 
Generically, this would not invalidate the previous argument made for a single Wilson- 
Dirac fermion and one would expect to have to flne tune. Strictly speaking, this is true for 
any flnite k, but in practice, any flne tuning (and the accompanying chirality violating 0{a) 
effects) could be ignored if the appropriate numerical coefficients were sufficiently small. 
Generically, one expects the coefficients to be of order one. The exponential suppression 
of these coefficients is therefore potentially very useful, but a bit mystifying. Some subtle 
cancelation must occur when one chooses the action in the way indicated by the overlap 
or, equivalently, either by the see-saw or by the F-N mechanisms. 

This paper has three main parts: In the first (next section) we show that in gauge 
backgrounds satisfying a certain criterion there is exponential mass suppression. We also 
show that instantons evade this suppression by violating the criterion. Thus we have 
a separation of gauge backgrounds into two classes, with approximate conservation of 
chiral symmetry in one and order one breaking in the other. The first class contains 



some neighborhood of the trivial gauge orbit. This structure mirrors our understanding in 
the continuum. However, this analysis does not shed light on the mysterious cancelation 
required in Feynman diagrams containing internal gauge field lines. Nor is it obvious to us 
that, even at one loop order, the gauge fields that violate the criterion occur indeed with 
zero probability. In section 3 we complete the picture by analyzing all the contributions 
to the zero momentum fermion propagator, (in Feynman gauge) to one loop order in 
perturbation theory. The analysis reveals the cancelations between different fiavors that 
are needed for the suppression to hold, indicates that other sets of orthogonal polynomials 
would also work just as they would at tree level and shows, as expected, that highly 
nonlocal pure gauge interactions could spoil the suppression. 

Mainly for future use we present in section 4 the structure of the free propagator in 
a way that should generalize to other sets of orthogonal polynomials. The full expression 
for the Chebyshev case (truncated overlap) is also given. Some of the technical details are 
relegated to an appendix. Our derivation is different from the one adopted in [6] (and later 
generalized to the truncated case in [8,12]) and leads to a particular form of presenting the 
result that is both more concise and, we feel, more insightful. We make no direct use of 
these results here, but they should be useful in other situations. 

The final section contains a summary and outlook. 



2. Light quark propagator in various gauge backgrounds. 



First, we establish our notation. It is very similar to that of [3]. The action S 
(appearing with a minus sign in the exponent when used to define the partition function) 
is: 

S=,:^Sn + SF + Sr,f. (2.1) 



1 
2~g 



T^^Sg + S'f + Spf. 



Sg is a pure gauge action whose detailed form is irrelevant here. The other two terms 

contain the fermions as Grassmann variables {Sp) and the pseudo-fermions as normal 

numbers, (Spf). We use both for the fermions and pseudo-fermions the following left-right 

structure: 

/ $1 \ 

$2 



, $2fc 
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X^ 

\xV 



(2.2) 



Xj ' are left or right Weyl fermions in the notation of [14]. Similarly one defines $s. 
Our convention is that vector gauge interaction appear diagonal in D. We suppressed all 
space-time, spinorial and gauge indices, displaying explicitly only the left-right character 
and flavor. 

The lattice is taken to have L^ sites. Our basic building blocks in the lattice Dirac 
matrix D will have size q x q where, in 4 dimensions, q = 2ncL'^- Uc is the dimension of 
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(2.3) 
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\ ... 

The matrix D is of size 2k x 2k and the entries are q x q blocks. 

The matrices B and C are dependent on the gauge background defined by the collec- 
tion of link matrices U^{x). The link matrices are of dimension Uc x Uc- /U labels the 4 
positive directions on a hypercubic lattice and U^ (x) is the unitary matrix associated with 
a link that points from the site x in the /i-direction. 
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{Bo)xaa,yf3b = ^<^a/3 ^[2(5xy(5a6 " 8y,x+jl{U^{x))ab — Sx,y+fi{UJ^{y))ab]- 



af{W^U,yb 



(2.4) 



M=l 



{B)xaa,yl3b — {Bo)xaa,yl3b + MQdxota,yl3b- 

The indices a, /J label spinor indices and take values 1 or 2. The indices a, b label color 
in the range 1 to Uc- The Euclidean 4x4 Dirac matrices 7^ are taken in the Weyl basis 
where their form is 



7m 











(2.5) 



As long as the parameter Mq satisfies < Mq < 1, S is positive for any gauge field 
background and the light tree level quark mass muL is small [3,5]: 
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(2.6) 
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A convenient choice for interpolating fields for the light quark ([6]) are the quantities 
xfiXfc (^1,^2) To obtain their correlation functions we introduce two source terms, X 



and y, in D [3 . 
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X and Y are qxq matrices which act as arbitrary sources for the bilinears Xiciai^)XkBbiy) ^ 
Xkaai^)xf(3biy) • '^o obtain the correlation functions we need to expand the ratio ^^etDp'/ 
around X = 0, y = to Unear order in X and Y. Actually, D^f = D{X = 1, y = 1). In 
the appendix of [3] the following identity was established: 



det D{X, Y) 
det DPf 



det 



-X 





-T- 



-Y 



det [1 + T- 



(2.11) 



Here: 

T is positive definite for any gauge field and det T = 1 . Defining tfie fiermitian matrix H 
by 

T = e-^. (2.13) 

we see tfiat tr(H) = 0. 

Tfie spectra of H and T are gauge invariant. For a background U^{x) gauge equivalent 
to U^{x) = 1, H fias q positive eigenvalues and q negative ones. Around zero, H fias a finite 
gap. Let us write U^{x) = e*^"^^*^^) wfiere tfie A^{x) are fixed by some local gauge condition 
(Feynman gauge - we fiave no problems witli Gribov copies liere as we simply assume tfiat 
A^{x) satisfies tfie gauge condition to identify ttie orbit). We see tfiat tlie structure of tfie 
spectrum could be preserved in some neighborhood of (7 = if the probability distribution 
for A^{x) (at leading order this probability distribution is g independent) suppresses large 
A^{x) strongly enough. Link configurations causing H to have equal numbers of positive 
and negative eigenstates, and with the minimal absolute value eigenvalue Cmin satisfying 
Gmink >> 1, are referred to as satisfying the "perturbative" criterion for k. For any k one 
can find a small enough range in g that a typical link will satisfy the perturbative criterion 
with high probability but it is difficult to make this statement more precise. This difficulty 
is the main reason for also performing a direct one loop calculation in the next section. 

To compute the two determinants in eq. (2.11) we use two different orthonormal 
bases, and represent the matrices by their matrix elements between the two bases. This 
trick produces expressions that are equivalent in form to those of the overlap ([1]), but 
now with the truncation effects made explicit. The first basis is denoted generically by v 
and is gauge field independent, v is indexed by a pair of indices: the first is either r or / 
while the second is J = 1,2, ....q. Each J is a short-hand for a triplet of indices: a Weyl 
index, a lattice site and a gauge group index (J = (a, x, a)). 

v^'^'^ = ("o") ' ""^'''^ ^ (^°0 ' ""^^^ ^ ^'^P^'^^^^y (2-14) 

The other basis, generically denoted by w, is made out of ortho-normalized eigenvectors 
of H. These vectors do depend on the gauge fields and transform covariantly under gauge 
transformations. Without restricting generality we can take the gauge group as SU{Nc), 
and therefore, once the bases are used for evaluating determinants, the gauge dependence 
disappears on account of the unimodularity of the group elements and the ultra-local 
structure of the basis v. 



We assume that the background satisfies our perturbative criterion. The q positive 
eigenvectors of H are denoted by w^'^' \ J = 1,2, ...q and the q negative ones by w^"^'"^': 



The 2q x 2q unitary matrix relating the v and the w bases is given by: 
U = 



(2.15) 






jjpr jjpl 
Tjnr TJ"^^ 



(2.16) 



We choose the phases so that det U^'^ = (det [/"^)* (this is possible since the unitarity of 
U implies Idett/^''! = jdett/""'!, as shown in [1]). With this choice detU = 1. For any 
2q X 2q matrix Z it is true then that det Z = det(^/)^(^j) w^^^^'^ Zv^'^'^^ where ^ = p,n and 
7] = l,r. 

Applying this to our basic expressions, we obtain: 
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On the other hand, we have, for large k, 



detfl + T-'^] -det 



.fcE? 



e^'^Hij 



5ij 

where we now used only the w basis on both sides. At large k we end up with 
det D{X, Y) 



(2.19) 



det 



-yj{P,I)\y{-r,J) yj{P,I)tYv^^'-^^ 

1X,i-riJ)'\ Xy{r,J) yj{-n,I)-\ y{l,J) 



det DPf 

Setting X = y = we obtain the overlap result [1]: 
det D{X = 0,Y = 0) 



(2.20) 



detL>(X = l,y = 1) 



det[/Pn^ = |det[/"^P. 



(2.21) 
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Just like in [3] , our derivation made no use of fermion creation and annihilation op- 
erators. The latter are essential for chiral gauge theories, and, as a consequence, can be 
also used in the vector-like case [1]. But operators techniques become quite clumsy when 
dealing with the truncation of the overlap, and unnecessarily complicate the derivations. 

Now, we need to expand to linear order in X and Y. But, if we replace any of X or y 
by zero, the dependence on the other variable drops out. We therefore conclude that, up 
to corrections exponentially smaU in k, < Xiaaix)Xkt3biy) >u and < Xkaai^)xff3biy) >u 
vanish. The same derivation would show that, if the number of positive and negative 
eigenvalues differs by two say, the above correlation functions do not have to vanish. Note 
that it is only one part of the perturbative criterion that gets violated. As long as the 
second part of the criterion is satisfied, namely an appropriately gapped spectrum, the 
large k limit, be it zero or nonzero, is approached exponentially fast. One could say that a 
weaker criterion, namely one which contains only the gap requirement, defines the subset 
of configurations to which semi-classical considerations apply [3] . 

As already mentioned, it is difficult to make precise probabilistic estimates on how 
often we have at most e^i^k ~ 1 for any given k and g. Therefore, the result of a 
perturbative calculation of the light quark propagator in Feynman gauge, even at one loop, 
is not predicted with certainty by the result of the present section. Had we been able to 
claim that gauge fields (assumed a priori to produce i7's having equal numbers of positive 
and negative eigenvalues), violating the e^i^k » 1 for a given fc, occur with vanishing 
probability as g is smaller than some small, ^-independent, e > 0, we could conclude, 
just on the basis of this section, that exponential mass suppression holds to all orders in 
perturbation theory. In the absence of this claim, the main importance of the present result 
is in establishing an effective decoupling between the left and right components of the light 
quark in "perturbative" backgrounds while also showing specifically how this decoupling 
gets spoiled in other backgrounds. 



3. Exponential mass suppression to one loop. 

The infinite mass matrix Ai was introduced in [6]. As emphasized there, although one 
wrote down a theory that formally looked vector-like, due to the impossibility of rotating 
all the fields so as to make Ai hermitian, the theory ended up being chiral. Of course, this 
slightly paradoxical situation occurred only because of the infinite dimensionality of Ai. 
With infinite mass matrices inherently non-hermitian, the free fermion propagator had to 
be separated into left and right parts, and separate expressions had to be written for each. 
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This approach has been taken over to the vector-hke case in [11] and [12]. Clearly, in the 
vector-like case it constitutes an unnecessary complication, particularly in the truncated 
overlap situation, where it should be transparently clear that all we are dealing with are 
several lattice Dirac fermions mixed in a special way. 

So, following [3], we change bases in flavor space to make the mass matrices hermitian. 
At the beginning, to avoid confusion, we shall denote explicitly the direct product structure 
between Dirac spinor space and the rest. Writing out the explicit index dependence on 
flavors {i,j = l,2...k) we have, in the new flavor basis [3], 



Dij = 7^ 



W„, 5.., + 1 



M "^iJ 



M, 



ij, 



(3.1) 



where the W^ were defined in (2.4) and a sum over // is implied. The hermitian matrix M 
has the following flavor structure: 



/ 




M 



\B 
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-1 
-1 B 













0/ 



(3.2) 



The entries are ricL'^ x ricL'^ matrices. 

To set up perturbation theory, set U^ (x) = 1 and go to Fourier space. The site depen- 
dence of M can be diagonalized and each entry becomes block diagonal with momentum 
{p) dependent ric x ric blocks. There are L^ such blocks in each entry. Actually, the factor 
in group space is unity, so we simply have each block of M represented by one real function 
of momentum. Let the matrices made out of these representatives, now of dimension kxk, 
be denoted by m{p). 

For any p the matrix rn?{p) is tridiagonal and has been diagonalized in [3]. Let us 
summarize what we need here. The eigenvalues are denoted by /^^(p) and are given by 



//2(p) = 1 + h''{p) - 2b{p)XM, Kp) = 5^(1 - cosp^) + Mo, 



(3.3) 



where the A^ (p) are the k real roots of the polynomial equation 

UkiX) = bip)Uk+iiX). 



(3.4) 



The Uj{X) are the Chebyshev orthogonal polynomials of the second kind and play a central 
role in what follows. We are quite convinced that other mass matrices, associated with 



11 



other sets of orthogonal polynomials, would work similarly and in order to stress that 
we shall try to avoid using any specific properties of the Chebyshev polynomials. The 
dependence on the momentum p only enters through b{p), but this could easily be changed 
in other variants. 

The eigenvectors corresponding to the above eigenvalues are proportional to Uj{Xs{p)): 



J2m\p)^JU,{Xs{p)) = filipmXsip)). 



(3.5) 



The orthonormal eigenvectors are given by 



where, for any A, 



0,{Xs{p)) = N{Xs{p))U,{Xs{p)) = N,{p)UjiXs{p)), 



N\X) 



E-=if//(A) 



(3.6) 



(3.7) 



For definiteness, we choose to define N{X) as the positive square root of (3.7). 

The orthogonal character of the polynomials ensures that all the A^ (p) are distinct for 
any fixed p. We can arrange the Ag (p) in descending order, so that s = 1 corresponds to the 
smallest mass. For each s, Xs{p) is a smooth function of p, as no crossings can occur. The 
associated set of k eigenvectors are globally and smoothly defined over momentum space. 
The lack of degeneracy for any p implies that diagonalizing 'w?{p) has also diagonalized 
m{p). Therefore, 



^mij{p)Oj{Xs{p)) = Hs{p)Oi{Xs{p)). 

The above equation defines the sign of ps (p) ■ More explicitly it can be read off from 

b{p) 1 



(3.8) 



f^s{p) 



Uk{Xs{p)) Uk+iiXs{p))' 



(3.9) 



or from 



f^s{p) 



1 



Hp) 



+ b{p)-2Xs{p) 



Uk{Xs{p)). 



(3.10) 



Ps{p) has the same sign as Uk{Xs{p))- From the above two relations one also derives 

-\ 2 



l + b^{p)-2b{p)Xs{p) 



Hp) 



UkiXsip)) 



(3.11) 



in agreement with (3.3). 
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It was shown in [3] that, as long as 

0<6(p)<l + p (3.12) 

Xi{p) is larger than unity. Let the region ([6]) of momentum space where < b{p) < 1 + -^ 
be denoted by R. For p E Rwe have, as fc ^ oo ([3]), and up to exponentially suppressed 
corrections, 



Mp) ~ 2 



bip) + 



Hp)\ ' 



f,t{p)r^b''{p)[l-b\p)r. (3.13) 



For all 2 < s < fc with any p and also for s = 1 with p ^ R \Xs{p)\ < 1 and Hs{p) 
stays finite and nonzero. Thus, the masses of all the heavy (doublers and extra fiavors) 
16/c — 1 Dirac fermions are given by the ultra-violet cutoff times order one coefficients. 
The corresponding eigenvector components, Oj(As(p)), are of typical order -7= and have, 
in general, an oscillatory behavior as a function of j. The normalization constants Ns{p) 
are also of order -7= in these cases. 

However, for s = 1 and p E R , where the light quark state resides, the components 
Oj(Ai(p)) do not oscillate but rather vary exponentially for large j (we assume k is large). 
Similarly, Ni{p) is exponentially large in that region. All this is a direct consequence 
of the exponential growth of the polynomials Uj{X) outside their interval of orthogonality 
(|A| < 1) and their boundedness within. This property will be shared by other sequences of 
orthogonal polynomials, at least as long as the interval of orthogonality is a finite segment. 

The above information is basically all we wish to use when analyzing the light fermion 
propagator to one loop. In momentum space the full free propagator is given by: 

Note that the matrices in the numerator and the denominator commute so there is no 
ordering problem and (3.14) is unambiguous. We ignored the trivial color dependence. 
Introducing the diagonal kx k matrices /x(p), with /ii(p), H2{p), ■■■IJ^kip) along the diagonal 
and the orthogonal matrices 0{p) made out of the eigenvectors Oj{\s{p)) as columns, we 
can write the free propagator as 

G^'\p) = 0(p)^^|^54/^0^(p). (3.15) 

p^ + /i^(p) 

Above, we dropped the explicit direct products and anything that is unity in the appro- 
priate space. 
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We now proceed with the calculation of the fermion propagator to one loop. As 
mentioned in the introduction we only wish to set things up and establish exponential 
suppression of the numerical coefficient of the 0{a~^) radiatively induced mass for s = 1 
and p E R, where a is the lattice spacing. We shall avoid any explicit calculation that 
is tangential to our goal. We choose to define the quark mass from the expansion of 
the inverse propagator around zero four momentum. Since all fermions are massive this 
expansion is not infrared divergent. 

Expanding (3.1) to order (/^ we have: 

Aj = Df^ + Mlf + g-i^,Wl^^5,, + gM^^^P,, + g^-i^W^^^S,, + g'^M^^^P,,. (3.16) 

Here, Pij = Si^k+i-j implements physical parity exchanging the left and right components 
of the light fermion. 

The propagator, G = -^, is also expanded to order g'^, and after that averaged over the 
gauge fields with a Gaussian measure, using Feynman gauge. The order g term averages 
to zero and we can write: 

where the self energy S is given by 

S = (I)W-^dW-Z)(^)), (3.18) 

with (...) denoting gauge averaging. In (3.18) we used the following short-hand notations: 

Dff=7,Wl,'^S,,+M('^P,, 



(3.19) 



Dff^jM'^5,j+M^'^P, 



We are interested in the self energy in momentum space, near zero momentum. Since 
we wish to find out the mass shift of the lightest mode (labeled by s = 1), we only need 
the diagonal matrix element of the self energy in the unperturbed light fiavor eigenstate, 
whose wave function in fiavor space is given by Oj{Xi{p)). 

First we wish to establish that we are not really interested in the vicinity of zero 
momentum, since all we would get from it is the wave function renormalization constant 
Z. By showing that Z cannot be exponentially large in k, so it cannot affect our conclusion 
about the possible exponential suppression of the light quark radiatively-induced mass, we 
can restrict our analysis to the self energy strictly at zero momentum. That Z cannot 
be exponentially large in k is quite obvious: The matrix O in (3.15) is orthogonal so its 
entries are bounded. There are no divergences worse than logarithmic in the infinite k 
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limit (when one of the quarks becomes massless). Thus, at most, Z could have a linear 
dependence on fc, and we can forget about Z altogether. 

The {D^ ">) ("tadpole") term has no free propagators. The factor containing W}i has 
no left-right terms (terms commuting with 75) so does not contribute. To estimate the 
contribution of M^^' we only need to compute the expectation value of P in the s = 1 
state. The structure of M^"^' is irrelevant, since all it determines is a multiplicative factor 
of order one. 

The following identity follows directly from the recursion relations for the Chebyshev 
polynomials: 

E U,(\)U,,,-,(P) = I'lil^zIhiM, (3,20) 

It implies, in particular, 

k 

J2 t/.(A)[/fc+i_,(A) = ^[/^+i(A), (3.21) 

i=l 

where prime denotes differentiation with respect to the argument. 

For the normalization Ns{X) we need the following identities, also directly derivable 
from the recursion relations for the orthogonal polynomials: 

fa.(A)£/.(p)4^^i±ii^l^«4^i^i±iM^^. (3.22) 

^-^^ 2 A — p 

This implies, in particular, 

1 2 

^'^^^ ^ Elif/f(A) ^ t/^+i(A)f/.(A)-[/.+i(A)[/^(A)- ^^-^^^ 

For the M'^ "> contribution we need: 

k 

^O,(Ai(0))Ofc+i_.(Ai(0)) = ^iVf(0)[/^+i(Ai(0)). (3.24) 

i=l 

For any A > 1 and large fc, Uk{\) ~ (2A)''~^ and A^(A) ~ (2\)k~i ■ Since Ai(0) > 1 we 
have proved the desired exponential suppression for this contribution. 

The first quantity in eq. (3.18) generates 8 contributions from the two terms in D^^' 
and the two terms in G^^' . Only four of these contributions commute with 75 and are of 
interest: 
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The VF(^)-VF(^) term will need the quantity defined below, 

k 

fwwip) = i2^I(p)-rMr-V (3.25) 

where 

k 

Ys{p) = Y,Oj{\imOAUp))- (3.26) 

The W^^'-M^^' and M^^'-W^^' terms are equal to each other and for them we need the 

quantity 

k ^ 

fwM{p) = fMwip) = y2Xs{p)Y,{p)^^ — J--, (3.27) 

where 

k 

Xs{p) = Y,OMimOk+i-j{\s{p)). (3.28) 

j=i 



For the M^^^-M^^) we need 

k 



s=l 



l^s{p) 

p^ + nl{p)' 



fMMip) =J2^'iP) -/_:2f y (3-29) 



Each of the f{p) functions goes into a loop integral over the four momentum p. Under 
the integral f{p) is multiplied by the explicit momentum dependence coming from the 
respective vertices and by the gauge propagator which we simply take as 

^ - ^ ■ p^ = 2sm^. (3.30) 



2E^(l-cosp^) p2' -M 2 

Overall factors coming from the Casimir in group space, 27r's in the integrals, etc., are 
irrelevant to us here. 

It is now that we are in a position to see why something "miraculous" has to happen. 
If exponential suppression is maintained, we expect every one of the contribution to be 
suppressed individually, since, for example, we can imagine varying M^^' without changing 
VF*^^), while maintaining the tree level mass hierarchy. Consider, for example, fww{p)- 
Y^(p) is bounded by 1 via the Schwartz inequality. For arbitrary p we would expect the 
bound to be an order of magnitude estimate, and then we would guess fwwip) to be of 
order k. The smallness of /ii(p) for p E R hardly seems to make a difi^erence. To get a 
suppression the alternating signs of the masses /Xs [p] must play a significant role, inducing 
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exponential suppression by almost perfect cancelations. Similar considerations apply to 
the two other terms. 

Until now we have made no use what-so-ever of the structure of the vertices. We 
expect to need to use the fact that M^^' is linear in p for small p, since it is related to the 
RG irrelevance of the Wilson mass term. 

Let us start with the VF(^)-VF(^'' term. From equations (3.4), (3.9), (3.22) we derive 
a simpler formula for Ns{p): 



NUp) _ 2 

f,s{p) b{p)U',^,{Xs{p)) - U',{Xs{p)) 



(3.31) 



Using now also equation (3.22) we obtain 



Ys'ip) 



2^s(p) 



iVi(0)^fc(Ai( 0)) 
6(0) 



n 2 



HO)-b{p) 
Ai(0)-A,(p) 



1 2 



b{p)u'{XM) - uj^iXsip)) 



(3.32) 
It is now evident that for s 7^ 1, Y^{p) is indeed order 1, as estimated above. The sign 
carried by the /J^sip) is also carried by the denominator of the last term in (3.32). Inserting 
into eq. (3.25), and using (3.3), we get, for p 7^ 0, 



fwwip) 



dz 



iVi(O) 
2 [^i{0) 



[b{p) 



1 



(3.33) 



c27Ti [Ai(0)-z]2 p^ + l + bHp)-2b{p)z b{p)Uk+i{z) - Uk{z) 



Here C encloses tightly all the zeros of (f){z,k,p) = b{p)Uk+i{z) — Uk{z). Along the real 
axis, (j){z,k,p) changes sign k times and the cancelations this causes are captured by 
deforming the contour of integration. Deforming the contour to a circle at infinity the line 
integral drops out and we pick up contributions from the two extra poles zi{k) = Ai(0) and 
Z2{p) = ol-bh) "'" Wl "'" Kp)]- Both zi{k) and Z2{p) are positive and larger than unity, so 



2lb{p) ^ b{p) 



JVi(O) 
^^1(0) 



is of 



when plugged into cf){z, k,p) give exponential suppression, as the prefactor, 

order one. Moreover, < zi{k) — 1 is bounded away from zero for all k, while < Z2{p) — 1 

is bounded away from zero for all p. 

However, the exponential suppression is wiped out at very small p'^ since, by the 
definition of Ai(0), (p{zi{0),k,0) = 0. The region of very small p where we have no 
exponential suppression is best analyzed using the original expressions (3.25) and (3.26). 
Since Ys{0) = Sis by ortho-normality, the dominating term at very small momenta is 
-2^ 2I -) ■ When multiplied by the propagator (3.30) we see that the suppression factor 
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//i(0) will be multiplied by log/U^(0). This behavior is expected from the continuum, and 
indeed the small p region of the loop integral is the place where continuum perturbation 
theory is reproduced ([15]). We learn that we should expect extra multiplicative factors 
of order k in front of the exponentially suppressed term. Here, we basically ignore these 
factors. The region of small p^ can be taken to extend up to p^ ~ |^i(0)|" with < 
a < 1 and the contribution from the rest of the integral is then seen to be exponentially 
suppressed. 

We now turn to the W^^^^-M^^^ and M^^^-W^-^^ terms. Manipulations similar to the 
above produce: 

^ ( .y( . NUO) b{p) - 5(0) Uk+i{Xim - Uk+i{Xs{p)) . . 

'^^^ '^^^ 2M0) MO) - XsipW b{p)u',^,iXM) - KiUp))' ^' ^ 

Again, we represent fwM{p) by 

fwM{p) = ^^^[b{p)-bm 

dz Uk+i{Xim-Uk+i{z) 1 1 



(3.35) 



Tc 2m [Ai(0) -z]^ p^ + l + 62(p) - 2b{p)z b{p)Uk+i{z) - Uk{z) ' 

The contour is first defined as above, and then deformed as before to infinity. Again 
the line integral does not contribute, and we end up only with contributions from the poles 
at Zi{k) and 2^2 (p)- From z = zi{k) and momenta not too small we obtain a contribution 
of order one from the integral and the prefactor ^l-J provides exponential suppression. 
From z = Z2{p) we obtain a contribution which, as a result of t/fc+i(Ai(0)) << Uk+i{z2{p)) 
for not too small p and k large enough, is also of order one. Thus the overall exponential 
suppression is, using (3.31), by 

NUO) 



m(o) 



[2zi{k)]-\ (3.36) 



For p ^ (3.27) will be dominated by the s = 1 term, given by Xi(0) -2,^i( ) ■ In 
addition we have a factor of J^ from the gauge field propagator and a factor going as p^ from 
the M^^' vertex combined with the 7^p^ term from G^^'{p). Thus, the region of momenta 
where (3.35) is not suppressed is too small to eliminate the exponential suppression found 
for momenta away from zero. Note that ^i(O) is also exponentially small. 

The last term to be analyzed is of type M^^>-M^^' . We need the following expression, 
derived with the help of (3.20), 



X^p) ^\Nnp)Nno) 



Uk+i{Xi{0))-Uk+i{K{p)) 



Ai(0) - Xsip) 



X!{0)=\[N!{0)UU,{X^{0))f. 
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(3.37) 



With it we obtain: 



fMM{p) = lNf{0)J2 



' Uk+i{\im-Uk+i{\s{p)) ]^ 

Ai(0)-A,(p) 
From (3.29) we see that fMM{p) has a finite hniit of order one as p^ -^ 0. In the Feynman 
diagram the gauge propagator singularity at p^ -^ is canceled by the momentum depen- 
dence coming from the M'^ "> vertices. Therefore, we can assume that p^ is sufficiently far 
from zero to permit individual treatment of the terms obtained from an expansion of the 
binomial squared in the numerator of the last factor in (3.38). There are three terms. The 
most "dangerous" contains the s-independent large constant t/|_|_]^(Ai(0)) (see below for 
the other two terms). The contribution from this term is: 

1 l + h''{p)-2h{p)z 1 ^^-^^^ 

[Ai(0) - z]2 p2 + 1 + h^{p) - 2b{p)z b{p)Uk+i{z) - Uk{z) ■ 
The prefactor is of order unity but the integral is exponentially suppressed as in the cases 
above. The other two terms, are also exponentially suppressed: The cross term has an 
exponentially small prefactor and an integral of order one. In the last term one needs 
to replace in the sum over s the factor t/|_|_]^(As(p)) by j^, ^2( )-2b{ )\ ( ) using equation 
(3.11). This avoids too fast growth at complex infinity in the associated contour integral. 
Now, one sees the exponential suppression easily. 

This concludes the argument. Note that a gauge propagator which is more singular 
in the infrared would have invalidated our argument by enhancing the relevance of the 
very low momentum region in the diagram. It would be hard to see this directly from the 
analysis of section 2, and this provides one example where the calculations of section 3 are 
seen to be explicitly needed. If one needs explicit, complete formulae the above analysis 
applies because all sums over s can be evaluated in closed form with the help of the complex 
contour integrals shown. A similar technique was used in [16]. All we do applies directly 
also to the "almost" supersymmetric case discussed in [3] by methods similar to [17]. 



4. Free fermion propagator. 

Using the methods of this paper we derive first a formula for the free propagator. We 
shall make the formula completely explicit for the Chebyshev case, but we wish to present 
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it in a less explicit manner first because, at that stage, the structure would be the same 
for other sets of orthogonal polynomials. 

According to (3.14) (see also [6]) we only need to invert the second order Dirac oper- 
ator. So we need the matrix 

^ 0(p)^^^^3— 0^(p). (4.1) 



Explicitly, 

A,(rt^i:iV|(rt ^-'«;y^'' . ,4.2) 

Although not indicated explicitly, A^^ also depends on k. Using the same techniques as 
before we arrive at 

A-(v) = 2b(v) I — ^-(^^^^-(^^ I (A 3) 

''^^^ ^^' fc 2m U,iz)[bip)U,+^{z) - U^iz)] p^ + 1 + b^p) - 2b{p)z- ^ ' ' 

Here the contour C encircles precisely only the roots of b{p)Uk+i{z) — Uk{z). Deforming 
to infinity we pick up one contribution from the pole at z = Z2 (p) and k — 1 contributions 
from the roots of Uk{z), denoted by xf- 



^ U',{xt)Uk+i{xt) p2 + 1 + b^{p) - 2b{p)xt ^^_^^ 

U,{Z2{P))U,{Z2{P)) 



Uu{z2{p)mp)Uk + l{z2{p)) - Uk{z2{p))\ 



The first term will be written as A9(p) while the second is, up to normalization, a pro- 
jection matrix on a state that, for p G -R, has entries exponentially decreasing with the 
distance of i from k. At k, the components are order unity. The second term represents 
the contribution of a state dominated by flavors near k. This state clearly is representing 
the almost massless quark. On the other hand A? {p) vanishes if either i or j is equal to 
k. Thus, it represents some other set of fc — 1 states. Jointly, these states span a /c — 1 
subspace in flavor space. This subspace is not dependent on the momentum p. 

Using the recursion relations we see that we can replace in the denominator the quan- 
tity Uk+i{xt) by —Uk-i{xt). Introducing now the orthogonal matrix 0^'^~^^ {p indepen- 
dent) defined as 



{k-i).2 1 (4.5) 
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we see that, for i,j = 1, ...k — 1, 
















fc-l ^{k-l)^{k-l) 

A^-(v)=\ '' '' = 


( ^ 


1 
u 










\p^ + ml{p), 




ml{p)=l + h'{p)-2h{p)3, 


J = 


/ 
1 



1 


1 




1 




. 
. 

1 . 


.. 
.. 
.. 

.. 

.. 1 


o\ 




1 

0/ 



(4.6) 



Above, the matrix J has dimensions (/c — l)x(/c — 1). The form of the J term in 
the matrix rMo(p) was determined by considering what matrix would have the eigenvalues 
Xt, t = l,2,...,/c— 1, by analogy with [3] where the inverse problem was solved. Note 
that although we defined the matrix TOq(p) it does not have a simple (sparse) square root, 
unlike m^ (p) . All k — 1 eigenvalues of TOq (p) are positive and bounded away from zero for 
all momenta p because all the zeros of Uk{z) are in the interval (—1, 1) for any k. 

It is clear that the above separation of the second order operator differentiates between 
the contributions of the light Wilson fermion and the other k — 1 heavy Wilson fermions. In 
fiavor space the heavy Wilson fermions span a fixed k — 1 dimensional subspace. The light 
Wilson fermion is associated with a direction that is almost orthogonal to the heavy sector 
of flavor space for p E R. The heavy doublers contained in the light Wilson fermion appear 
to have unsuppressed overlaps with the heavy Wilson fermions, but some cancelations 
prevent them from giving substantial mass to the light component of the light Wilson 
fermion. 

The full propagator is easily obtainable now, once the explicit form of m{p) is used. 

To this point we only did manipulations that we believe would generalize to other 
sets of orthonormal polynomials. We now write down the answers allowing ourselves to 
use specific properties of the Chebyshev polynomials. In [12] similar quantities where 
calculated following steps laid out in [8] who generalized the method of [6] to the finite k 
case. 

Following [6] we introduce the positive quantity a{p): 



Z2{P) 



P 



+ 



1 



+ h{p) 



= cosh(a(p)), h{p) 



Mo + ^p' 



< Mo < 1. (4.7) 



_h{p) b{p) 

It is important to note ([6]) that a{p) is smooth and bounded away from zero for all 
momenta p. Using standard manipulations of trigonometric identities and contour inte- 
grations, we derive in the appendix the following explicit expression for the free propagator. 
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including all finite k effects: 

sinli{a(p)[niin(A; 



^^Jip) 



hk- j)]} sinh{a(p)[min(i, j)]} 



b{p) sinh[a(p)] sinh[/ca(p)] 
sinh[zQ;(p)] sinh[jtt(p)] 



(4.8) 



smh[ka{p)]{b{p) smh[{k + l)a{p)] — smh[ka{p)]} 
The term on the second line in (4.8) contains the light particle. One large k limit is 
obtained by taking k to infinity with k — i = i' and k — j = j' held fixed. In this limit, 
i',j' = 0, l,2....oo, and we obtain, for p^ ^ 0, 



A 



(l)oo 



*J 



(p) 



*(p)K'-i'l 



-a(p)(i'+j') 



»(p)(j'+j') 



+ 



Note that 6(0) exp[a(0)] = 1, but when written as above, the pole at pp' - 



(4.9) 



evident. As /c ^ oo, for p^ = Q and fixed, finite i', /, the second term diverges as ' ^ 



becomes 

2fc 



Ko) 
refiecting the exponential smallness of the mass. Another large k limit is obtained taking 

/c ^ oo with i and j kept finite, i,j = 1,2, The second term in (4.8) disappears for 

p^ 7^ while the first is similar to (4.9) because of its invariance under i ^ k — i, j — ^ k—j: 

g-a(p)|i-j| _ g-a(p)(i+j) 



(2)c 



*J 



Again, for p^ = 0, special analysis is required. In the two above limits we choose to keep 
two distinct groups of particles. These two groups decouple at infinite k. Exponential 
mass suppression requires this decoupling to hold approximatively also at finite k. In both 
cases the limits p^ ^ and A; ^ oo do not commute. The regime where this lack of 
commutativity is felt can eventually become important if the gauge propagator is replaced 
by something with a stronger divergence in the infrared. For the ordinary singularity, in 
four dimensions, the lack of commutativity appears to have no major effect to one loop 
order. However, in two dimensions, the ordinary singularity is already sufficiently strong 
to require a deeper investigation. 

The full Dirac propagator is given by: 

(7(0) (p) = A{p)[-i-f^p^ + m{p)] = [-i-f^p^ + m{p)]A{p), (4.11) 

where the k x k matrix m{p) has the following structure: 



/ 




m 












-1 



-1 b{p) 
\b{p) 



-1 

b{p) 








b{p)\ 






/ 



(4.12) 
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If one wishes to keep a finite mass at /c = oo it lias to be introduced by hand [1,8]. 
Since it is obvious that the flavor k strongly mixes with the light fermion, one can use 
this component of the fermions as an interpolating field for the light particle [6] . Even at 
finite k one may find some advantages to keep such an explicit extra mass term [3] . In our 
notation it amounts to replacing the mass matrix m{p) by to(p) + /iN, where n is the new 
mass parameter and N has a single nonzero entry: N^j = SikSjk, so as to couple only the 
left and right components of the interpolating field. 

In equation (4.9) we have the propagator for ^u, = 0. It is easy to obtain the propagator 
for/x^O, 4°^(p): 

^1^^P^^ + m{p) + ^N 
Since N obeys 

N2 = N, NG(°)(p)N = £?o(p)N, (4.14) 

where go{p) = Gkl{p), one can simply expand in /i and resum the series: 

G'(f)(p) = G^''\p)-^lG^''\p)^G^''\p)+^?G^''\p)go{p)^G^''\p)- 

^^G^^\p)gl{p)^G^^{p) + = G'W(p) - GW(p)--^^NG(°)(p). ^^'^^^ 

Note that as far as spinorial indices go, go{p) is still a matrix, and ordering is important. 
So, the free propagator with the extra mass term p is given by: 

gT ip). = G^ ip) - G^ ip) ^MTT^S^ (^) • (4-16) 



*j 



Since t/jk is the interpolating field for the light fermion one can imagine integrating out all 
the other fermion fields [6]. The action, for /U = 0, would obviously be tpk—jw^ipk- Now the 
introduction of /x 7^ has a trivial effect, merely adding the term nt(jkil^k- Therefore, the 
new free propagator for the V'/c field is: 

GT(p)kk = ^^^- (4-17) 

It is easy to see that this coincides with eq. (4.16). G),l{p) has a relatively simple form: 

^(0) . ^ ^ -^1^^P^^ s\nh{ka{p)] + h{p) sinh[a(p)] 
^^ ^^' b{p) smh[{k + l)a{p)] - smh[ka{p)] ' ^ ' ' 

It is unclear at the present whether employment of the free propagator given above 
would result in a simplified analysis in the truncated overlap case. Our hope is that the 
more general treatment we presented in section 3 will promote a search for more efficient 
mass suppression schemes. Note that the analysis in section 2 relied on the particular 
structure of the truncated overlap from an earlier stage. Nevertheless, it too can be made 
more general, using the formulae given in the appendix of [3]. 
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5. Summary and Outlook. 

The evidence in favor of exponential suppression of the light quark mass in the trun- 
cated overlap and potentially in related systems is fairly strong. Thus there is a new 
competitor to existing methods designed to reduce/eliminate 0{a) effects in QCD sim- 
ulations. In practice, the requirement Mq > cannot be met at practical values of the 
gauge coupling. It remains to be understood whether the required negative values of Mq, 
and the potentially associated "exceptional configurations" are better, worse or similar to 
those affecting ordinary Wilson fermions. In [3] it was suggested that for — 1 < Mq < 
and odd /c's one would be dealing with the 9 = it regime of QCD. It would be interesting 
to see how much of the present paper extends to this case. Whether this suggestion should 
extend to the k = 1 case (ordinary Wilson fermions) or not is too early to even guess. 
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6. Appendix. 

This appendix deals exclusively with the Chebyshev case. The results therefore are 
not expected to generalize easily to other sets of orthogonal polynomials. Of course, the 
truncated overlap (domain wall) is covered. 

For Chebyshev polynomials we have: 

Uk{x) = '-^^^; x = cos(^). (Al) 

Hence, 

xt = cos{et); et = ty- t = l,....,A;-l. (^.2) 

k 

A short calculation then gives: 



Ni =--r——- -— = - sm {9t), 

J U^{xt)Uk-i[xt) k 

Of^Of,-'^ =1 [cos((z - m) - cos((^ + j)et)] . 



{A.3) 



Thus, 

Q 1 ^ COS (^t7r) — cos (^tyr) 

''^^' ^ kbip) f~[ 2 cos [if] - 2cosh(a(p)) ' ^ '^' 
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Replace t hj k — t in the above sum and average the two expressions. 
For i + j even we obtain 



1 



A«,(p)=2cosh(a(p))-— ^5^ 



cos (^y^tTr) — cos (^^t7r) 



kbip) j-[ [2 cos (f )]' - [2cosh(a(p))]^ 



(A5) 



Note that we have extended the sum to t = k because the corresponding term in the sum 
vanishes. Introduce now two no n- negative integers, /, /': 



/ 



i+j 



I' 



\^-J\ 



2 ' 2 

A'^j{p) is seen to be expressible in terms of a quantity gi{p), 

k _ _ _ /r, , / \ / k 



{A.Q) 



9i 



(p)-E 



cos 



(2^4) 



J-i 



t [2cos(f)] -[2cosh(ab))]^ 



Real(^- 



[zt - e2«(p)] [zt - e-2«(p)] J ' 

{A.7) 
where zt = e^™'i^ , t = 1, ...,k are all the solutions of z^ = 1. 

Consider now a complex line integral / over a closed circle at infinity: 



dz kz 



k-i 



1 



1 



_ e2"(p) 



27ri z^ — 1 
I vanishes for any / > 0. This leads to 

k cosh[{k-2l)a{p)\ 



z - e-2°(p) 



{A.l 



9i{p) 



2 sinh[2a(p)] sinh[/ca(p)] 

from which A^Ap) can be obtained for even i + j. 
For odd i + j we define 



{A.9) 



I 



i+j + 1 



I' 



K-j| + i 



Note that now 1,1' > 1. A few lines lead us to 

1 



Kip) 



kb{p) 



[gi{p) + gi-i{p) - gv{p) - gv-i{p)\. 



(A.IO) 



(All) 



Collecting all results we see that the expressions for even i -\- j and odd i + j are the 
same and our final result becomes: 



A?,(p) 



cosh [(A; — \i — i|)o(p)] — cosh [{k — i — i)a{p)] 



2b{p) sinh[«(p)] sinh[/c«(p)] 
A few more steps produce equation (4.8). 



(^.12) 
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